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Abstract 

Anisotropies of Young's modulus E, the shear modulus G, and Pois- 
son's ratio v of all 2D symmetry systems are studied. The shear modulus 
and Poisson's ratio of 2D crystals have fourfold symmetry. Simple neces- 
sary and sufficient conditions on their elastic compliances are derived to 
identify if any of these crystals is completely auxetic, non-auxetic or aux- 
etic. Examples of all types of auxetic properties of crystals of oblique and 
rectangular symmetry are presented. Particular attention is paid to 2D 
crystals of quadratic symmetry. All mechanically stable quadratic crys- 
tals are characterized by three parameters belonging to a prism with the 
stability triangle in the base. Regions in the stability triangle in which 
quadratic materials are completely auxetic, non-auxetic, and auxetic are 
established. 

1 Introduction 

Mechanical properties of elastic media are described by the bulk, Young's, and 
the shear modules. Poisson's ratio is another important characteristics of such 
materials. With the exception of bulk modulus, generally, these characteristics 
are anisotropic, i.e. they depend on direction n of load and direction m of lateral 
strain. The Poisson ratio is especially interesting because it can be negative. 
Materials which exhibit a negative Poisson's ratio are referred to as auxetic 
media. Auxetics respond to imposed uniaxial tension with lateral extension in 
place of expected contraction. Such media can find interesting applications in 
future technologies. Ting and Barnett introduced classification of the auxetic 
behavior of anisotropic linear media [1] . 

In our previous papers we derived expressions for the inverse of Young's 
modulus E(n), the inverse of the shear modulus G(m, n) and Poisson's ratio 
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^(m, n) of all 2D [2] and of 3D crystalline media of high and medium symmetry 
[3]. In another paper [4] we studied the above mechanical properties of all 
mechanically stable cubic media. To achieve this goal, we used a set parameters 
introduced by Every [5]. All stable cubic elastic materials are represented by 
points of a prism with the stability triangle in the base. We established regions 
in the stability triangle in which 3D cubic elastic materials, in accord to Ting 
and Barnett classification [1], are completely auxetic, non-auxetic, and auxetic. 

The aim of the present paper is to advance the study of anisotropy properties 
of the mentioned material constants of 2D crystals and formulate necessary and 
sufficient conditions on elastic compliances to identify whether any given 2D 
material is completely auxetic, auxetic or non-auxetic. The quadratic materials 
will be studied in detail. We will show that their elastic properties can be 
studied in the same manner as those of 3D cubic materials. For this purpose 
we shall use parameters, which were 2D analogues of parameters introduced by 
Every [5] (cf. Appendix A). 

To explain microscopic mechanisms underlying the property of auxeticity, 
various models are proposed. Two-dimensional microscopic crystalline models 
make an interesting class of auxetics (cf. [6]- [10]). Each model of such kind 
belongs to one of 2D crystal systems, therefore its auxetic properties are char- 
acterized by Poisson's ratio suitable for this system. The classic example of 
such structure is re-entrant honeycomb structure shown in Fig. 1 (cf. rcf. [7]). 
In Sect. 4.1 we show examples of crystals of oblique and rectangular symmetry 
which are completely auxetic, auxetic and non-auxetic. 



honeycomb-rect . png 



Figure 1: Left panel: Re-entrant honeycomb structure, right panel: the corre- 
sponding 2D lattice is rectangular centered. 



2 Anisotropy of material constants of 2D media 

The treatment of anisotropy of Young's and shear modules and of Poisson's ratio 
is based on expressions derived in the companion paper [2]. These expressions 
do not depend on the choice of the reference axes. However, here we prefer to use 



2 



their form valid in crystalline reference axes (CRA). We denote the unit vectors 
parallel to the crystal axes by u, and v. We shall outline the derivation of 
explicit angular dependence of the elastic material constants for oblique crystal 
system. Suitable expressions for rectangular, quadratic and isotropic media 
can be obtained by reduction the number of independent and non-vanishing 
compliances. 

In the case of oblique 2D materials we obtained [2] 

B _1 (n) = s n (nu) 4 + s 22 (nv) 4 + 2 (s 33 + s 12 ) (nu) 2 (nv) 2 

+2V2s 13 (nu) 3 (nv) + 2\/2s 23 (nu)(nv) 3 , (1) 
-u(m, n)/E(n) = sn(mu) 2 (nu) 2 + s 22 (mv) 2 (nv) 2 
+2s3 3 (mu)(mv)(nu)(nv) + s i2 [(mu) 2 (nv) 2 + (mv) 2 (nu) 2 ] 
+V2si3 [(mu) 2 (nu)(nv) + (nu) 2 (mu)(mv)] 
+V2s 23 [(mv) 2 (nu)(nv) + (nv) 2 (mu)(mv)] , (2) 

[4G(m, n)]" 1 = Sll (mu) 2 (nu) 2 + s 22 (mv) 2 (nv) 2 
+s 33 [(mu)(nv) + (mv)(nu)] 2 /2 + 2si 2 (mu)(nu)(mv)(nv) 
+V2 [si 3 (mu)(nu) + s 23 (mv)(nv)] [(mu)(nv) + (mv)(nu)] , (3) 

with Sij — 1, 2, 3) is given in the companion paper [2]. 

The vectors of stretch n and strain m are two mutually perpendicular vectors 
of the unit length. Their components fulfill three conditions 

(nu) 2 + (nv) 2 = 1, (mu) 2 + (mv) 2 = 1, (nu)(mu) + (nv)(mv) = 0. (4) 

This means that only one component of n and m is independent. In CRA the 
angle ip between the vector u and vector n is defined by the scalar product 
(nu) = cos <p. The components of m obey the relations (mu) = mi = — n 2 = 
— sin ip, (mv) = m 2 = n\ — cos ip. 

2.1 Anisotropy of Poisson's ratio 

In CRA the function [—i/ (m, n)/E (n)] can be written in the following form 

-v (ip)/E (ip) = [-v o /E o ] + r v cos {Aip + ip v ) , (5) 

where 

[-v o /E o ] = [(S n + S 22 - S e6 ) /2 + 3,S 12 ] /4, 

r v = sJ{S 2& - S 16 f + [S 12 - (S n + S 22 - S 66 ) /2] 2 /4, 

^26 — Sie . . 



For rectangular, quadratic and isotropic media S\e = S 2 6 = 0. This means 
that for them ip v = 0. Adapting Eq. (5) for oblique materials to rectangular, 
quadratic and isotropic materials we obtain 



+ [S12 - (Sn + S 22 - S 66 ) /2] cos(4^)/4, (7) 



E T ((f ) \ E T J 
where 

(-v T /E T ) = [(S n + S 22 - S 66 ) /2 + 3S 12 ] /4. (8) 
For quadratic media Sn = S22 and 

-v^/E^ip) = (-is q /E q ) Q + [S12 - (2Sn - S 66 ) /2] cos(4^)/4, (9) 

where 

(- V£q) = [(25n - S 66 ) /2 + 3S 12 ] /4. (10) 
For isotropic media Sqq = 2 (Sn — S12) (cf. [2]), hence form Eq. (9) we obtain 

u is /E is = -S 12 = G 12 /(G n - C? 2 ). (11) 

From last of Eqs. (17) and Eq. (11) it follows that v is = — Si 2 /Sn. An isotropic 
material is mechanically stable if Sn > |Si 2 | (cf. [2]). Thus, in agreement with 
rcf. [11], two-dimensional isotropic systems can exhibit Poisson's ratio within 
the range —1 < z/, s < +1. 

2.2 Anisotropy of the shear modulus 

For oblique materials in the CRA Eq. (3) yields the following expression for 
(4G)-\ 

[4G (¥>)] _1 = (4G o ) 1 + r G cos(4<^ + Vg), (12) 

where 

(4Go) _1 = (Sn + S 22 - 2S 12 + S 66 )/8, 
(S 66 + 2S 12 - S n - S 22 ) 2 /4 + (S 26 - S 16 ) 2 1 V ' /4, 



r G = 



2 (Si 6 — S 2 6) 

tan^c = 



(Sae + 2S12 — Sn — S22) 
For rectangular, quadratic and isotropic media, i/jq = and 

[4G r (^)]" 1 = (4GO)" 1 + (S 66 + 2S 12 - S n - S 22 ) cos(4^)/8, (13) 

and for quadratic media 

^Gq^)]" 1 = (4G q )" 1 + (S 66 + 2S 12 - 2S n ) cos(4^)/8, (14) 

where 

(4G q )- 1 = (2S n - 2S 12 + S 66 ) /8. (15) 
For isotropic media, the last formula gives the familiar result [2] 

Gr/=2/(C 11 -C 12 )>0. (16) 
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2.3 Anisotropy of Young's modulus 

Young's modulus has more complicated structure than Poisson's ratio and the 
shear modulus. The expressions given in Sect. 2 yield 

E^ 1 = Sn cos 4 p + Sn sin 4 p + 2S\e cos 3 p sin p + 2Sie cos p sin 3 <p 

+ (25*12 + S 6e ) cos 2 ip sin 2 tp, 
E~ x = Sn cos 4 p + 5 22 sin 4 p + {2S 12 + 5 66 ) cos 2 p sin 2 p, 
E^ 1 = Su (cos 4 p + sin 4 + (25 i2 + See) sin 2 p cos 2 <p, 

Er^Sn- (17) 

Formulas (7)-(17) can be obtained with the help of general formulas (1)- 
(3) for material constants, or directly in CRA from their definitions (cf. Eqs. 
(33)-(35) of ref. [2]). 

3 Auxetic properties of 2D media 

3.1 Oblique and rectangular crystals 

Since is positive, Eq. (5) implies that an oblique crystal is completely 
auxetic, i.e. v is negative for all direction pairs (m,n), or for all values of the 
angle p, if, and only, if 

{-"o/Eo)o > 6 o) > 0. (18) 

If 

{-Vo/E o ) < 0, | {-v o /E o ) | > r(°), (19) 

a crystal is non-auxtetic, i.e. ^(n, m) is positive for all direction pairs (m, n) or 
for all angles p. 
If 

| (-v /E )o I < r v , (20) 
a crystal is auxetic, i.e. there can be found auxetic direction pairs (m,n) (angles 
(p) for which the right hand side of Eq. (5) is positive, and non-auxetic direction 
pairs for which it is negative. 

Unlike in the case of 2D crystals of oblique symmetry, for crystals of rect- 
angular symmetry, the parameter can be negative. Anisotropy of elastic 
material constants and auxetic properties of quadratic crystals will be studied 
in the next section. 

3.2 Auxetic properties of quadratic crystals 

In the case of quadratic crystals, one may use eigenvalues cj, cl, and cm in place 
of three elastic constants Cn, C12, and Cee, but the more reasonable choice is 
the use of three parameters, which are the analogs of the parameters introduced 
by Every [5], namely, Si = (Cn + Cee) an d two dimensionless parameters 

S2 = {C 11 -C 6e )/s 1 , s 3 = K/ Sl , (21) 
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where K — (Cn — Ci 2 — Cqq) /2 is the elastic anisotropy parameter. One may 
express the elastic constants Cu (I, J = 1,2,6) by s\, s 2 and s 3 parameters, 
namely 



di = Sl (1 + s 2 ) /2, C12 - si (3s 2 /2 - s 3 - 1/2) , C 66 = *i (1 - s 2 ) • (22) 



These parameters parametrize phase velocities of quadratic crystals in a natural 
way (cf. Appendix A). 

In terms of parameters si s 2 , and s 3 a quadratic material is stable, if s\ > 
and s 2 , S3 belong to the stability triangle (ST for short) in the plane (s 2 , S3) 
(cf. Fig. 2) defined by three inequalities 



Sides of ST are defined by lines AB (s 3 = 2s 2 ), BD (s 2 = 1) and AD (s 3 = 
s 2 -l). 



triangle2d.png 



Figure 2: Triangle ABD of stability (ST) of 2D crystals of quadratic symmetry. 
Points (si, s 2 , S3) characterizing mechanically stable quadratic materials belong 
to the prism (si > 0), the base of which is the stability triangle. Points of the 
interval < s 2 < 1 represent isotropic media. 

The set of inequalities 



1 > 82, (1 - S 2 + S3) > 0, (82 - 8 3 /2) > 0. 



(23) 



(-!/,/£;,)„> 0, (-u g /E q ) >\r^\ 



(24) 



G 



No. 


Compliances 


Every's parameters 


region of ST 


I 


S12 > 


|«2 - «3 - ^ < 


above AC 


II 


Su + S 12 - S 66 /2 > 


s 3 - 3s 2 + 1 > 


above BF 


III 


ri q) > 


S 3 > 


upper part of ST 


IV 


S 12 < 


|s 2 - s 3 - i > 


below AC 


V 


Su + S 12 - S 66 /2 < 


s 3 - 3s 2 + 1 < 


below BF 


VI 


ri q) < 


S3 < 
Table 1: 


lower part of ST 



is equivalent to inequalities shown in the I and II row of the Table 1, and defines 
the region of ST in form of the triangle AEB representing crystals of quadratic 
symmetry which are completely auxetic. 

A crystal of quadratic symmetry is non-auxetic when 

(-u g /E q ) <0, \ -(v q /E q ) Q \> \r v \. (25) 

The inequalities (25) are equivalent to inequalities IV and V of the Table 1 , and 
circumscribe the quadrilateral region FECD of ST. 

A crystal of quadratic symmetry is auxetic, provided that 

\{-to q /E g ) Q \ < |r(?>|. (26) 
If > 0, then the inequality (26) is equivalent to two pairs of conditions 

HV£,) D >0,r('»>(-4^)o. 
{-to q /E q ) Q < 0, r(?> > - (-^ q /E q ) , (27) 

which are equivalent to inequalities I and V defining the triangle EBC . 
If rl°^ < 0, then in agreement with inequalities I and IV of Table 1 

(-4v q /E q ) >0,(-4v q /E q ) <-ri*\ 
{-Av q /E q ) < 0, - (-4v q /E q ) < -r(?>, (28) 

bound the triangle AFE. 

We shall mention that points of the interval S3 = 0, < s 2 < 1/3 represent 
isotropic completely auxetic 2D materials, whereas points of the interval S3 = 
0, 1/3 < s 2 < 1 represents isotropic non-auxetic 2D elastic media. 



4 Anisotropy of elastic material constants of 2D 
crystals 

4.1 Anisotropy of elastic material constants of crystals of 
oblique and rectangular symmetry 

In the case of crystals of oblique and rectangular symmetry we restrict the 
discussion of the anisotropy to Poisson's coefficient. Using compliances collected 
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in the Appendix B we draw polar plots for Poisson's coefficient of oblique and 
rectangular media. They are shown in Figs. 3 and 4. 



poisson_obl_cz j as . png 



Figure 3: Dependence of ^(m, n) /E(n) of crystals of oblique symmetry on angle 
ip for materials with compliances collected in Table 2. (a) and (b) - auxetic, (c) 
- completely auxetic, and (d) - non-auxetic behavior. 

Inspecting these figures we conclude that there exits completely auxetic, 
auxetic and non-auxetic crystals of oblique and rectangular symmetry. 

4.2 Anisotropy of elastic material constants of crystals of 
quadratic symmetry 

We shall characterize the response of the elastic quadratic material to uniaxial 
tension in a direction n (nn = 1) by the dimensionless Young's modulus e(n) = 



where S'^ kl = s\Sij k i- Repetition of a suffix in a product of tensors implies the 
usual summation with the respect to that suffix over the values 1,2. In Fig. 5 
we present the dependence of dimensionless Young modulus e = E/s\ on the 
parameters S2, S3 for the same characteristic directions as above. 

Similarly, the response to the shear stress will be characterized by the di- 
mensionless shear modulus 17(11, m) = G(n, m)/si, (mm = 1, nm = 0) 



E(n)/s! (cf. ref. [4] and [2]) 



e 1 (n) = (n ® n) • S' • (n ® n) = riinjS' ijkl n k n h 



(29) 
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(m ® n) ■ S" ■ (m <g> n) = m i njS' iokl 'm k ni. 



(30) 



4ff(m, n) 
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poisson_rect_czjas .png 



Figure 4: Dependence of ^(m, n)/_E(n) of crystals of rectangular symmetry on 
the angle ip for materials listed in Table 3. (a) - non-auxetic, (b) - completely 
auxetic, (c) and (d) - auxetic behavior. 
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young . png 



Figure 5: Quadratic crystals. Contour map of dimensionless Young modulus 

e q (s 2 , s 3 ; n) for 

n =< 1, 1 > and n =< 1, > 
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The dependence of g q on S2 and S3 for the same characteristic directions as 
in the case of Young's modulus is shown in Fig. 6. 



shear . png 



Figure 6: Quadratic crystals. Contour map of dimensionless shear modulus 
9 q (s2, S3; m, n) for n =< 1, 1 >, m =< — 1, 1 >, and n =< 1, >, m =< 
0,1 > 

The dimensionless Poisson's ratio z/(n, m) is the ratio of the lateral contrac- 
tion along direction m to the longitudinal extension along direction n perpen- 
dicular to m 

i/(m,n) = -e(n)(m®m) ■ S' ■ (n ® n) = ' ljkl (31) 

Tit ?lu £>tuvw TLw 

For quadratic crystals the dependence of dimensionless Poisson's ratio v/s\ 
on parameters S2, S3 for two characteristic directions n =< 1,1 > and n =< 
1,0 > is shown in Fig. 7. Inspecting these figures we note that in the case 
of quadratic materials v obeys the inequality established for isotropic media in 
Sect. 2.1. 

Finally, the results of our calculation of angular dependence of Poisson's ratio 
for points belonging to the characteristic regions of ST are presented in Fig. 8. 
The solid line represents complete auxetics (the region AEBG in Fig. 2). The 
dotted line represents complete non-auxetics (the region FECD) and the broken 
line - non-complete auxetics (the region AFE and ECB). For them Poisson's 
ratio has to be anisotropic. The point S2 = 0.1, S3 = represents complete 
isotropic auxetics. All these contours have fourfold symmetry in agreement 
with formula (9), but their orientations depend on coordinates S2, S3. 
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poisson.png 



Figure 7: Quadratic crystals. Contour map of dimensionless Poisson's ratio 
v q (s2, S3; m, n) for n =< 1, 1 >, m =< — 1, 1 >, and n =< 1, >, n =< 0, 1 > 

A Phase velocities for quadratic materials 

Suppose that an acoustic wave propagates in a quadratic medium with the 
density p in the direction n. In the reference frame related to the symmetry 
axes of a quadratic lattice, elements = n p Ci P j q n q of Christoffel's matrix T 
are given by [2] 



Ci, C2 , and the anisotropy parameter K are the analogues of parameters intro- 
duced by Every [5] 



C\ — (C11 + Cqq) , C2 = (Cn — Cqq) , K = (C11 — C12 — 2C66) • (33) 



The eigenvalues of T are equal to pvj (J = ±1), where j enumerates two 
phase velocities. Denote by T the trace of Christoffel's matrix. For crystals of 
quadratic symmetry T = TrT = (Cn +C^). Upon making the replacement 
2pv 2 = S + T [5], one arrives at the following eigenvalue problem equation for S 



Tn = (dm? + C'eent) = C 2 n\ + C 66 , 
r 22 = {Ciml + C 66 n?) = C 2 n2 + C 66 , 
T 2 i = (C12 + C 66 ) nin 2 = (C 2 - if) nin 2 . 



(32) 



|A-S/ 2 | 



2C 2 (n? - 1/2) -5 (C 2 - X) nm 2 
(C* 2 - if) nin 2 2C 2 {n\ - 1/2) - S 1 



= 0. 



(34) 
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Figure 8: Dependence of v q (m, r\)/E q (p) on the angle ip for several points be- 
longing to characteristic regions of ST for quadratic crystals. Left diagram - 
full line represents auxetic behavior, broken line - non-auxetic. Right diagram 
illustrates completely auxetic behavior of anisotropic (full line) and isotropic 
material (broken line). 
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On 


<Sl2 


S22 


Sie 


S26 


5*66 


aux. propert. 


(a) 


0.105 


0.024 


0.372 


-0.071 


0.311 


0.626 


auxetic 


(b) 


0.037 


-0.006 


0.633 


0.038 


-0.015 


0.090 


auxetic 


(c) 


0.178 


0.111 


0.823 


-0.057 


-0.053 


0.601 


compl. auxetic 


(d) 


0.037 


-0.001 


0.321 


0.033 


0.0004 


0.910 


non-auxetic 



Table 2: Compliances for crystals of oblique symmetry. They represent materi- 
als with all types auxetic properties. 





Su 


5*12 


5*22 


£66 


aux. propert. 


(a) 


0.077 


-0.008 


0.215 


0.798 


non-auxctic 


(b) 


0.100 


0.007 


0.315 


0.277 


compl. auxetic 


(c) 


0.030 


0.019 


0.672 


0.928 


auxetic 


(d) 


0.100 


-0.080 


0.971 


0.368 


auxetic 



Table 3: Compliances for crystals of rectangular symmetry. They represent 
materials with all types of auxetic properties. 



On expanding the determinant in Eq. (34) one obtains a quadratic equation, 
whose solutions reads 



1 

Yp 



d ± 2\j {C 2 - Kf n\n\ - Cf (n\ - 1/2) (n\ - 1/2) 



(35) 



For isotropic media, Cqq = (Cn — C12) /2 and See = 2 (Su — £12), hence 
for them K = 0. Thus, the condition for elastic isotropy is that K = 0. When 
this is satisfied, one obtains the phase velocities for 2D isotropic media [2] 



v+ = VI = y/Cu/p, V- EE v t = V (Cu - C12) /2p. 



(36) 



B Examples of compliances for crystals of oblique| 
and rectangular symmetry 

We have written a computer program which finds numerical values of compli- 
ances obeying the mechanical stability conditions described in ref. [4]. Several 
such sets for oblique and rectangular materials are collected in tables 2 and 3 

These values are used for calculating angular dependence of Poisson's coef- 
ficient for oblique and rectangular materials (cf. Sect. 4.1). 



C Instabilities on the borders of ST 

Inspecting Figs. 5 and 6 one notice that for considered directions n and m 
Young's modulus and the shear modulus G(m, n) vanish on various borders of 
the stability triangle (cf. Table 4). The bulk modulus n is defined as [14] 

n- 1 = h ■ S ■ I 2 - S Ujj = S n + 2S 12 + S 2 2- (37) 
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Figure 9: Dependence of dimensionless bulk modulus k' on S2 and S3 for 
quadratic crystals, k' vanishes on the upper side of ST 
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E 


n = (1,1) 


: BF 


n = (1,0); AB, AD 


G 


n = (1,1), m = 


-1,1); AD 


n = (1,0), m = (0,1); BD 



Table 4: Instabilities of modules for various directions. Sides of ST on which 
instabilities occur are indicated. 

For quadratic crystals k = [2 (S n + S^)}' 1 = (Cn + C 12 ) /2 = cj/2 = si (2s 2 - s 3 ) /2 =| 
si/t' vanishes on upper side of ST (cf. also Fig. 9). Vanishing of modules marks 
phase transitions. The vanishing the shear modulus is accompanied by soft 
modes. For 3D cubic crystals we discussed all such instabilities and possible 
low symmetry phases in ref. [15]. For such crystals instability related to n = 
leads to the discontinuous isomorphic phase transition. 
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